Thermal Desorption of Hydrogen Prom Graphane 



L. A. Openov* and A. I. Podlivaev 

Moscow Engineering Physics Institute (State University), 115409 Moscow, Russia 
* E-mail: LAOpenovCimephi.ru 

ABSTRACT 

The process of hydrogen desorption from graphane (graphene sheet saturated by hydro- 
gen adsorbed from both sides) has been studied using the method of molecular dynamics. 
The temperature dependences of the time of desorption onset for various hydrogen cover- 
ages on graphene are calculated and the corresponding activation energies in the Arrhenius 
equation are determined. It is established that graphane exhibits a rather high thermal 
stability that makes possible its use in two-dimensional electronics even at room temper- 
ature. For the same reason, graphane can hardly be considered as a promising hydrogen 
storage material for fuel cells. 
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Graphene, representing a monolayer sheet of carbon atoms, draws the attention of both 
basic science (as a medium featm^ing massless Dirac fermions) and practical applications 
(as a material for nanoelectronics) [1]. It is also of interest to study the properties of 
graphene derivatives (e.g., in the form of nanoribbons [2]), including those obtained by its 
chemical modification with polymers, nitrogen, hydrogen, etc. Sofo et al. [3] theoretically 
predicted the existence of graphane, representing a graphene sheet saturated by hydrogen 
adsorbed from both sides, and quite recently this prediction was experimentally confirmed 
[4]. 

Graphane (unlike graphene) is a dielectric and, in principle, it can also be used in na- 
noelectronics. The relatively large content (7.7 mass %) of hydrogen in graphane suggests 
that it probably can be used in hydrogen power engineering. In assessing the possibilities 
of the practical usage of graphane, it is necessary first to study its thermal stability with 
respect to hydrogen desorption and thus determine the range of permissible temperatures 
for potential applications. This Letter presents the results of an investigation of the des- 
orption of hydrogen from graphane, which was numerically simulated using the method 
of molecular dynamics. 

The initial sample of graphene was modeled by a hexagonally packed monolayer frag- 
ment consisting of 54 carbon atoms. In order to reduce the influence of the finite size of 
the sample (edge effects), the dangling bonds of sp- hybridized edge carbon atoms were 
saturated by hydrogen. This fragment was converted into graphane by attaching one 
hydrogen atoms to each carbon atom alternatively from the two sides of the monolayer. 
As a result, we obtained a C54H72 model cluster of graphane schematically depicted in 
Fig. 1. In this cluster, all atoms occur in the sp^-hybridized state [4] (in contrast to 
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the sp^-hybridization in graphene) and are shifted in the transverse direction (off-plane) 
by 0.2 - 0.3 A depending on the position relative to the center of the cluster. The C-H 
and H-H bond lengths in the graphane fragment are 1.10 and 1.50 - 1.53 A, respectively. 
Note that it is necessary to distinguish between 18 hydrogen atoms used to passivate the 
sample edge and 54 hydrogen atoms forming the graphane structure: the former H atoms 
are present in both graphene and graphane models, whereas the latter are present only 
in graphane; in macroscopic graphane samples of increasing size, the relative fraction of 
the former H atoms tends to zero, while that of the latter tends to unity. 

The aim of this study was to determine the temperature dependence of the time nec- 
essary for the detachment of one H atom (of the total of 54 hydrogen atoms entering into 
the graphane structure) from the C54H72 cluster. This corresponds to the removal of a 
1/54: 2 % fraction of H atoms from a macroscopic sample of graphane, that is, to the 
formation of a relatively large number of defects (sufficient to significantly change the 
electrical characteristics) . This task was solved using the method of molecular dynamics 
[5, 6]. At the initial moment, all atoms were randomly imparted velocities and displace- 
ments so that the total momentum and moment of the entire cluster were zero. Then, 
the forces acting on atoms were calculated and the classical Newton equations of motion 
were numerically integrated using the velocity Verlet method with a temporal step of 
to = 2.72 X 10"^^ s. The total system energy (a sum of the potential and kinetic energies) 
remained unchanged, which corresponded to a microcanonical ensemble (a system not in 
thermal equilibrium with the environment) [5, 6]. The relative motion of atoms was char- 
acterized by the dynamic temperature T that was determined from the following relation: 
(Ekin) = ^kBT{3n — 6) where {E^in) is the kinetic energy averaged over the time, fc^ is 
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the Boltzmann constant, and n is the number of atoms in the system. This approach was 
successfully used previously for the qualitative explanation of experimental data on the 
fragmentation of Ceo fuUerene [6] . 

An important question is how to calculate the forces acting on atoms at each step of 
molecular dynamics. In fact, this procedure reduces to calculating the total system en- 
ergy as a function of the coordinates of component atoms. Solving this task by the ab 
initio methods requires considerable computational facilities and allows the evolution of a 
system involving ~ 100 atoms to be traced for a rather short time on the order of ~ 10 ps, 
which is insufficient to accumulate a necessary statistics. For this reason, we employed a 
nonorthogonal tight- binding model [7], which reasonably compromise between the strict 
ab initio methods and oversimplified approaches based on classical interatomic interac- 
tion potentials. The proposed model rather well describes both small carbon clusters 
(e.g. fuUerenes [7]) and hydrocarbon structures (e.g., cubane CgHg [8]) and macroscopic 
systems [7]. We believe that the method of molecular dynamics in combination with a 
tight-binding potential provides a quite adequate description of the thermal stability of 
graphane, similar to the results obtained by applying the same approach to cubane [8]. 

Evolution of the C54H72 cluster was traced to the moment of detachment of the first H 
atom at various temperatures in the range of T = 1300 — 3000 K. The results presented in 
Fig. 2 show that as the temperature decreases, the desorption onset time r exponentially 
increases from ~ 0.1 ps to ~ 10 ns. The plot of ln(T) versus 1/T is well approximated by 
a straight line, which indicates that the process of hydrogen desorption can be described 
in terms of the following standard Arrhenius relationship: 



T-^{T) = A-exp 
4 



Ea 



ksT. 



(1) 



where A is the temperature-independent (or weakly dependent) frequency factor and £"„ 
is the activation energy (which can be determined from the slope of the straight line in 
Fig. 2). Statistical analysis of the results of modehng yields Ea — (2.46 ± 0.17) eV and 
A = (2.1 ±0.5) X 10^^ s~^. By substituting these Ea and A values into formula (1), we can 
evaluate r at low temperatures, which are inaccessible for computer simulations because 
of extremely long necessary computation times. At T — 300 K, this estimation yields a 
macroscopic value of r ~ 10^^ s, which is evidence for the possibility of using graphane in 
nanoelectronic devices operating at room temperature. As the temperature increases to 
600 K, the value of r drops to ~ 1000 s. These results do not contradict the recent data 
[4] , according to which the complete desorption of hydrogen from graphane was achieved 
by annealing it for one day in argon at T = 700 K. It should be noted that, in view 
of the exponential character of relation (1), any refinement of the Ea value can lead to 
rather significant changes in the values of r. Therefore, the above values of r at different 
temperatures should be considered as rough estimates. 

A detailed analysis of evolution of the model C54H72 cluster showed that the desorption 
of hydrogen from various parts of the cluster takes place at approximately equal proba- 
bihty . This fact agrees with the results of calculations of the heights U of energy barriers 
for the desorption of H atoms. These values were determined by studying the potential 
energy hypersurface of the given cluster as a function of the coordinates of component 
atoms (see [9]). The barrier height U amounts to 2.4 eV at the center of the cluster and 
increases to 2.7 eV at the periphery. As expected, these values are close to the desorption 
activation energies Ea- It should be noted that, in the course of simulations, we observed 
two cases of the desorption of H2 molecules rather than atomic hydrogen. 
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In order to study the dependence of the rate of desorption on the degree of graphene 
coverage by hydrogen (graphane corresponds to the total coverage, whereby the number 
Nff of nonpassivating hydrogen atoms is equal to the number Nq of carbon atoms in the 
sample) , we have also determined r as a function of T for two other systems representing 
graphene with a half coverage (^h — Nc/2, C54H45 cluster) and with a single adsorbed 
H atom (Nij — 1, C54H19 cluster). The corresponding values of the desorption activation 
energy {Ea = (1.74 ± 0.17) and Ea = (1.86 ± 0.14) eV, respectively) proved to be smaller 
compared to the value obtained for graphane, in agreement with the lower barriers for 
desorption {U = (1.2 — 1.9) eV, increasing toward the periphery) and with the general 
tendency of a decrease in the thermodynamic stability of hydrocarbons with decreasing 
relative content of hydrogen [3] . These results indicate that the time of the total desorption 
of hydrogen from graphane is determined by the initial stage, in which the hydrogen 
content decreases by (1 - 10) %. 

It is interesting to note that, in simulating the C54H45 and C54H19 clusters, we have 
repeatedly observed the fracture of the carbon framework much before the hydrogen 
desorption event. This behavior agrees with the data of experiments [4], in which an 
increase in the temperature of graphane annealing above a certain level led to damage of 
the graphene framework. Taking into account that (i) this effect is not observed in the 
C54H72 cluster and (ii) the rate of hydrogen desorption increases with decreasing content 
of H atoms, it may be recommended to perform the initial stage of annealing at a higher 
temperature and then slightly cool the sample. 

In conclusion, let us briefly consider the possibility of using graphane in fuel cells for 
electric car engines. In addition to meeting the requirements of high hydrogen content 
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(more than 6 mass %) and stability at room temperature, it is necessary that the fuel- 
cell material would be capable of rapidly (within ~ 1 s) and almost completely release 
hydrogen at a temperature not exceeding 400 K [10]. According to the results of our 
investigation, graphane does not obey the last requirement. This is probably related to 
the strong covalent C-H bonds which, on the one hand, ensure a high thermal stability 
of hydrogen chemisorbed on carbon nanostructures [10] and, on the other hand, sharply 
decrease the rate of hydrogen desorption. Therefore, the most promising directions of the 
practical use of graphane are apparently related to nanoelectronics. In this context, it will 
be of interest to study the phenomenon of thermoactivated migration of hydrogen through 
the graphene/graphane interface. Should this boundary remain atomically sharp at not 
very low temperatures, it will be possible to attempt at creating various nanoelectronic 
devices by selectively adsorbing hydrogen on graphene sheets and nanoribbons. 
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Fig. 1. C54H72 cluster modeling a graphane fragment. Large and small balls represent 
carbon and hydrogen atoms, respectively. 
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Fig. 2. Plot of the logarithm of desorption time r [s] of one H atom from the C54H72 
cluster versus inverse temperature T [K] . Points present the results of model calculations, 
solid line shows the linear approximation by least squares. 
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